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Impedance  measurements  of  both  electrodes  of  a  flooded  OEM  SLI  battery  at  various  SOC  and  with  various 
direct  currents  have  been  measured.  For  each  part  of  the  impedance  spectra,  an  electrochemical  process 
is  proposed  and  implemented  in  a  simulation  model.  By  simulation  of  impedance  spectroscopy,  spectra 
of  the  model  are  obtained  and  compared  with  the  measurement.  In  the  first  part,  the  focus  was  put  on 
inductive  semicircles,  while  in  this  second  part,  concentration  limitation  and  its  appearance  in  impedance 
spectra  is  investigated. 
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1.  Introduction 

In  the  first  part  of  this  paper,  inductive  semicircles  occurring 
in  the  spectra  of  both  positive  and  negative  electrode  of  lead- 
acid  batteries  have  been  regarded.  All  spectra  shown  in  that  part 
have  been  corrected  for  points  that  are  not  consistent  according  to 
Kramers-Kronig  consistency  check. 

Strictly  speaking,  the  impedance  can  only  be  calculated  if  three 
conditions  are  fulfiled:  linearity,  time  invariance  and  causality. 
Minor  violations  of  linearity  and  time  invariance  are  tolerable, 
which  is  then  called  quasi-linear  or  quasi-stationary.  They  can  be 
compensated  during  the  measurement  and  data  processing  [4]. 
If  the  violations  are  within  tolerable  limits  can  be  tested  with 
the  Hilbert  or  Kramers-Kronig  transform  [1,6,7].  Especially  non¬ 
linearity  is  not  a  problem  if  both  real  and  imaginary  part  changes  by 
the  same  factor  [7].  The  principle  of  Kramers-Kronig  is  to  calculate 
the  real  part  from  the  imaginary  part  or  vice  versa: 
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The  calculation  can  then  be  compared  with  the  measurement. 
From  Eqs.  (1)  and  (2)  it  becomes  clear  that  also  measurements  at 
co  =  0  and  co  oo  are  needed  to  calculate  the  integrals.  Schiller  et  al. 
[3]  have  developed  a  method  called  Z-HIT  to  overcome  this  prob¬ 
lem.  The  absolute  value  of  the  impedance  |Z|  is  estimated  from  the 
phase  angle  0  taking  into  account  the  frequency  range  from  starting 
frequency  cos  to  the  currently  regarded  frequency  co: 
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The  first  term  on  the  right  hand  side  is  the  logarithmic  Hilbert 
transform  and  the  second  term  is  a  correction  term  accounting  for 
the  limited  frequency  range. 

Now,  also  the  points  that  are  not  consistent  are  considered.  Sim¬ 
ulations  assuming  concentration  limitation  as  the  origin  for  the 
deformation  are  shown  and  compared  with  the  measurements. 


2.  Impedance  spectra  of  lead-acid  battery  electrode 

The  impedance  spectra  of  the  positive  and  negative  electrode 
shown  in  Figs.  1  and  2  [Part  I  of  this  paper]  are  corrected  for 
impedance  values  that  did  not  pass  the  Kramers-Kronig  consis¬ 
tency  check,  indicating  that  at  least  one  prerequisite  (linearity, 
stationarity,  causality)  was  not  fulfiled.  Figs.  3  and  4  show  the 
same  spectra  as  Figs.  1  and  2,  respectively,  but  including  those 
values.  Both  at  the  negative  and  the  positive  electrode,  the 
consistency  check  fails  at  low  frequencies  at  the  end  of  each 
measurement. 
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Fig.  1.  Impedance  spectra  during  charging  (left)  and  discharging  (right)  of  the  negative  electrode  of  a  60  Ah  flooded  lead-acid  battery  at  different  SOC  at  25  °C,  superposed 
direct  current  is  ±0.5  ho  and  the  SOC  was  adjusted  by  discharging.  Measured  frequencies  are  between  6  kHz  and  3  mHz  with  eight  frequencies  per  decade.  Points  that  do  not 
pass  Kramers-Kronig  consistency  check  are  not  shown;  the  complete  spectra  are  depicted  in  Fig.  3. 


Fig.  2.  Impedance  spectra  during  charging  (left)  and  discharging  (right)  of  the  positive  electrode  of  a  60  Ah  flooded  lead-acid  battery  at  different  SOC  at  25  °C,  superposed 
direct  current  is  ±0.5  ho-  The  SOC  was  adjusted  by  discharging.  Measured  frequencies  are  between  61<Hz  and  3  mHz  with  eight  frequencies  per  decade.  Points  that  do  not 
pass  Kramers-Kronig  consistency  check  are  not  shown;  the  complete  spectra  are  depicted  in  Fig.  4. 


Fig.  3.  Impedance  spectra  of  the  negative  electrode  during  charging  (left)  and  discharging  (right)  as  in  Fig.  1  (different  SOCs  at  25  °C,  superposed  direct  current  is  ±0.5  ho), 
but  including  impedance  values  that  do  not  pass  Kramers-Kronig  consistency  check.  Note  the  different  scaling.  Measured  frequencies  are  between  6  kHz  and  3  mHz  with 
eight  frequencies  per  decade  [2]. 
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Fig.  4.  Impedance  spectra  of  the  positive  electrode  during  charging  (left)  and  discharging  (right)  as  in  Fig.  2  (different  SOCs  at  25  °C,  superposed  direct  current  is  ±0.5  J2 o), 
but  including  impedance  values  that  do  not  pass  Kramers-Kronig  consistency  check.  Note  the  different  scaling.  Measured  frequencies  are  between  6  kHz  and  3  mHz  with 
eight  frequencies  per  decade  [2]. 
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Fig.  5.  Measured  and  calculated  (Z-HIT  reconstruction)  absolute  value  of  the  impedance  left:  negative  electrode  at  20%  SOC  during  discharge  with  0.5  ho  (corresponding  to 
Fig.  3  right),  right:  positive  electrode  at  20%  SOC  during  charge  with  0.5  ho  (corresponding  to  Fig.  4  left). 


Fig.  5  exemplarily  shows  a  comparison  between  the  measured 
and  reconstructed  absolute  value  of  the  impedance  for  selected 
spectra  of  the  positive  and  negative  electrode.  At  the  negative  elec¬ 
trode,  the  agreement  is  quite  good  down  to  a  frequency  of  about 
20  mHz  and  only  down  to  about  1 00  mHz  at  the  positive  electrode. 

A  failure  of  the  consistency  check  does  not  necessarily  mean  that 
the  measurement  equipment  does  not  work  well.  As  the  problem 
always  arises  at  low  frequencies,  it  is  more  likely  that  the  battery 
leaves  the  tolerable  range  of  deviation  from  linearity,  stationarity 
and  causality  after  some  time  of  constant  current  charging  or  dis¬ 
charging.  Consequently,  the  prerequisites  for  the  Kramers-Kronig 
transform  are  not  fulfiled  at  the  corresponding  frequencies. 

Most  probably,  the  concentration  limitation  of  the  reaction 
partners  Pb2+  and  S042-  ions  leads  to  instationary  and  acausal 
behaviour  of  the  battery.  Concentration  influences  charge  transfer 
in  a  way  that  small  concentration  leads  to  a  large  charge-transfer 
impedance.  At  the  negative  electrode,  concentration  limitation 
additionally  leads  to  “curls”  at  lower  states  of  charge  (Fig.  3  right) 
and  at  the  positive  electrode  at  90%  SOC  also  a  kind  of  curl  can  be 
observed  (Fig.  4  left)  and  the  real  part  of  the  impedance  becomes 
negative  in  some  cases.  Even  though  the  remaining  points  passed 
Kramers-Kronig  check,  this  does  not  mean  that  they  are  not  influ¬ 
enced  by  concentration  limitation.  Especially  the  values  close  to 
those  left  out  are  most  probably  also  influenced  by  concentration, 


but  the  variation  from  one  measurement  to  the  next  is  so  small  that 
the  measurement  can  be  assumed  to  be  quasi-stationary. 

3.  Simulation 

In  the  first  part  of  this  paper,  concentration  limitation  has  been 
neglected.  In  reality,  especially  at  low  frequencies,  the  limitation  of 
charge  carriers  can  be  observed  from  curls  or  large  impedance  val¬ 
ues  at  low  frequencies.  Instationarities  can  occur  at  both  electrodes 
and  both  during  charging  and  discharging,  which  means  that  a  lim¬ 
itation  of  both  Pb2+  ions  during  charging  (reverse  reaction  of  step 
2  in  Fig.  6)  and  of  S042-  ions  during  discharging  (step  3  in  Fig.  6). 
Such  behaviour  can  also  be  simulated  by  limiting  the  concentra¬ 
tion.  For  this,  frequency-dependent  factors  are  introduced  in  the 
differential  equations  of  the  negative  (Section  3.1 )  and  the  positive 
(Section  3.2)  electrode.1  Equations  to  describe  the  concentration  of 
Pb2+  ions  and  H2S04  have  been  derived  by  Sauer  and  simplified  by 


1  In  reality,  the  concentration  is  not  frequency  dependent.  However,  as  the  con¬ 

centration  falls  or  rises  with  proceeding  charge  or  discharge,  this  leads  to  an 

apparent  frequency  dependency  in  a  spectrum  with  superposed  direct  current.  For 

simplicity,  in  the  simulations  the  concentration  is  kept  constant  for  each  frequency, 
while  in  reality  it  changes  gradually  all  the  time. 
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Fig.  6.  Course  of  reaction  during  discharge  (see  Part  I  of  this  paper):  (1)  charge 
transfer,  (2)  desorption,  and  (3)  chemical  reaction. 

Thele  [5].  They  are  further  simplified  in  the  following  as  the  aim 
here  is  only  to  demonstrate  the  general  impact  of  concentration 
limitation. 

A  failure  in  Kramers-Kronig  consistency  check  (Z-HIT)  does  not 
necessarily  indicate  a  failure  of  the  measurement  device,  but  simply 
that  at  least  one  of  the  prerequisites  linearity,  causality  and  station¬ 
ary  is  violated.  The  simulations  shown  in  the  following  indicate 
that  concentration  limitation  leads  to  a  failure  in  the  consistency 
check,  so  it  can  be  assumed  that  this  is  also  what  happens  during  the 
measurements.  A  concentration  variation  with  frequency  violates 
at  least  the  demand  for  stationarity.  Stationarity  requires  that  no 
other  parameter  than  frequency  changes  from  one  measurement  to 
another.  Also  causality  is  not  fulfiled  because  the  impedance  is  not 
only  influenced  by  the  input  parameters  current  and  voltage,  but 
also  by  concentration,  which  cannot  be  influenced  from  outside. 

3.1.  Negative  electrode 

Concentration  factors  cH2so4(/)/cS2So4  and  cpb2+(f)/c£b2+  are 
added  in  the  adsorption/desorption  part  of  the  coverage  equation 
of  the  negative  electrode  (cp.  Eq.  (5)  in  Part  I  of  this  paper): 
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Fig.  8.  Kramers-Kronig  consistency  check  (Z-HIT)  for  the  simulation  with  4  l2 o  in 
Fig.  7. 


A  simplified  equation  to  calculate  the  H2S04  concentration  in 
the  pores  of  the  negative  electrode  is: 

CH2S04(t)  =  £ne^t)  J  SnegW  '  £>H2S04(Ch2S04) - ^2^  ' 

+  GMR,neg  +  GsR?neg^f  (5) 

£neg  is  the  porosity  of  the  negative  electrode,  ex  is  an  exponent, 
DH2so4  is  the  diffusion  parameter  of  sulphuric  acid,  Ax  is  the  diffu¬ 
sion  distance  and  GMr  and  GSr  are  generation/consumption  terms 
for  main  and  side  reactions.  During  impedance-spectroscopy  sim¬ 
ulation  it  is  assumed  that  the  SOC  and  thus  porosity  stays  more  or 
less  constant.  Gassing  is  not  considered  in  the  model  and  it  can  also 
be  assumed  that  the  diffusion  parameter  stays  constant.  Thus,  for 
a  rough  estimation  of  the  concentration,  assuming  further  that  the 
concentration  is  constant  during  the  simulation  of  each  frequency 
fi ,  Eq.  (5)  simplifies  to: 


^ f 1  =  feet  •  ((1  -  0(t))  •  exp  (jffv(t))  -  0(t)  ■  exp  R“)n%(t))) 
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Fig.  7.  Simulation  of  negative  electrode  with  frequency-dependent  concentration  limitation,  frequency  range  10  kHz  to  0.1  mHz,  eight  frequencies  per  decade,  with 
R0  =  0.8m£2,  boating  =  1  m£2,  Ccoating  =  100F,  a  =  0.6,  fi  =  0.5,  i'o  =  15A,  Cdl  =  1000F,  kct  =  0.001s_1  and  kdeSad  =  0.005  s-1,  left:  original,  right:  after  removal  of  points  that  fail 
Kramers-Kronig  consistency  check  (Z-HIT). 
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Fig.  9.  Simulations  of  positive  electrode  with  frequency-dependent  concentration  limitation,  frequency  range  1 0  kHz  to  1  mHz,  eight  frequencies  per  decade,  with  R0  =  0.1  m£2, 
Cdl  =  5000  F,  a  =  0.4,  J0  =  2  A;  left:  original,  right:  after  removal  of  points  that  fail  Kramers-Kronig  consistency  check  (Z-HIT). 


Ch2S04C fi)  =  CH2S04Cfi-l) 

|  £neg  '  Ph2SQ4  '  (CCh2SQ4(/i-1  )  ~  cfl2so4 )/  A*2 )  +  CMR.neg  ^ 

£neg 

With  GMr  =  1(2 10+  -  1  )/2F  and  dividing  by  the  equilibrium  con¬ 
centration  c^2So4  this  becomes: 


removing  the  values  that  did  not  pass  Kramers-Kronig  consistency 
check.  Apart  from  the  -1 12 o  simulation,  no  inductive  semicircle  is 
visible  any  more,  which  makes  it  practically  impossible  to  identify 
the  corresponding  parameters.  Comparison  with  the  corresponding 
simulation  without  concentration  limitation  [Part  I  of  this  paper] 
shows  that  the  remaining  values  are  still  partly  influenced  by  the 
limited  concentration. 
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For  i  =  1,  cH2so4(/i-i)/cg2So4  =  1  is  assumed. 

The  concentration  of  the  Pb2+  ions  can  be  calculated  from: 

dCPb2+W 

— ^ - =Q-s  (8) 

g  =  -//(2FVuni)  is  the  amount  of  ions  generated  or  consumed  by 
the  main  reaction;  the  negative  sign  is  needed  because  charging 
consumes  ions  and  discharging  generates  them.  Vuni  is  a  unified 
volume,  s  is  the  amount  of  Pb2+  ions  that  combines  with  SO42-  ions 
to  lead  sulphate.  It  depends  on  the  radius  distribution /(r,t)  of  the 
existing  sulphate  crystals: 

i  s /  r3  /(r’t)dr  (9) 


3.2.  Positive  electrode 

Neglecting  the  high-frequency  inductive  semicircle  and  assum¬ 
ing  that  no  adsorption/desorption  process  occurs,  the  faradaic 
impedance  of  the  positive  electrode  simply  consists  of  a  charge- 
transfer  resistance.  The  overall  impedance  is  then  a  purely  ohmic 
resistance  in  parallel  to  the  double  layer  capacitance  plus  an 
additional  series  inductance  L.  However,  the  Kramers-Kronig  con¬ 
sistency  check  of  the  measured  spectra  (Fig.  5)  suggests  that  strong 
concentration  limitation  is  present.  Similar  to  the  negative  elec¬ 
trode,  frequency-dependent  concentration  factors  can  be  added; 
here  to  the  Butler-Volmer  equation: 


Eq.  (9)  can  be  simplified  for  the  case  that  all  crystals  have  the 
same  radius.  This  simplification  is  good  enough  for  a  rough  estima¬ 
tion  of  the  concentration: 

s  =  4  •  n  •  N  •  r  •  Dpb 2+  •  (cpb2+  -  c°b2+)  (10) 

Assuming  again  that  the  concentration  is  constant  for  each  fre¬ 
quency  Eq.  (8)  becomes: 


(12) 
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Fig.  7  left  shows  simulated  spectra  with  concentration  lim¬ 
itation.  In  analogy  to  the  measurements,  high  frequencies  are 
simulated  first  so  that  limitation  becomes  only  visible  at  low  fre¬ 
quencies.  As  in  the  measurements,  curls  occur  at  low  frequencies. 
Fig.  8  shows  a  comparison  of  the  simulated  and  the  reconstructed 
absolute  value  of  the  impedance  corresponding  to  the  spectrum 
at  4/2q  in  Fig.  7  left.  Fig.  7  right  shows  the  simulated  spectra  after 


The  factors  are  calculated  analogously  to  Eqs.  (7)  and  (11),  but 
with  the  corresponding  values  for  the  positive  electrode  for  poros¬ 
ity  and  diffusion  distance  Ax.  The  generation  term  of  the  main 
reaction  Gmr  is  different  as  well: 


GMR,pos  — 


(3-2t+).J 


(13) 


2F 
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Fig.  1 0.  Kramers-Kronig  consistency  check  (Z-HIT)  for  simulation  with  4 12 o  in  Fig.  9. 

Fig.  9  left  shows  simulated  spectra  of  a  positive  electrode  with 
concentration  limitation.  Similar  to  the  measurements,  the  semi¬ 
circle  shows  a  “bump”  and  the  real  part  even  becomes  negative 
at  some  currents.  Fig.  9  right  shows  the  same  simulation  leaving 
out  the  points  that  do  not  pass  Kramers-Kronig  consistency  check 
(Fig.  10).  As  for  the  measurements  (Fig.  2),  not  much  of  the  spec¬ 
trum  is  left  and  it  is  hard  to  extract  impedance  parameters  from 
this.  Also,  it  can  be  seen  that  even  though  the  remaining  impedance 
values  pass  the  consistency  check,  they  are  still  affected  by  the  con¬ 
centration  limitation,  indicated  by  the  bend  toward  the  imaginary 
axis. 

4.  Conclusion 

Deformations  of  the  spectra  of  both  electrodes  of  lead-acid  bat¬ 
teries  due  to  concentration  limitation  have  been  discussed.  It  has 


been  demonstrated  that  a  frequency-dependent  concentration  of 
Pb2+  ions  and  H2SO4  leads  to  curls  in  the  spectra  of  the  negative 
electrode  at  low  frequencies  and  to  bumps  in  the  capacitive  semi¬ 
circle  of  the  positive  electrode  at  intermediate  frequencies,  similar 
to  the  measured  spectra.  Also,  Kramers-Kronig  transform  fails  for 
those  frequency  ranges,  just  as  for  the  measurements.  It  is  thus 
very  likely  that  concentration  limitation  is  the  cause  for  stationarity 
violations  and  consequent  deformations  of  the  spectra. 
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